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ABSTRACT 

The earliest generation of stars and black holes must have established an early 'Lyman- 
Werner' background (LWB) at high redshift, prior to the epoch of reionization. Because 
of the long mean free path of photons with energies hv < 13.6eV, the LWB was 
nearly uniform. However, some variation in the LWB is expected due to the discrete 
nature of the sources, and their highly clustered spatial distribution. In this paper, we 
compute the probability distribution function (PDF) of the LW flux that irradiates 
dark matter (DM) halos collapsing at high-redshift (z « 10). Our model accounts for 
(i) the clustering of DM halos, (ii) Poisson fluctuations in the number of corresponding 
star forming galaxies, and (iii) scatter in the LW luminosity produced by halos of a 
given mass (calibrated using local observations). We find that > 99% of the DM halos 
are illuminated by a LW flux within a factor of 2 of the global mean value. However, a 
small fraction, 10^* — 10~^, of DM halos with virial temperatures T^ir >W^ K have 
a close luminous neighbor within <10 kpc, and are exposed to a LW flux exceeding the 
global mean by a factor of > 20, or to J2i,lw > 10^ (in units of 10~^^ erg s~^ Hz~^ sr~^ 
cm~^). This large LW flux can photo-dissociate H2 molecules in the gas collapsing due 
to atomic cooling in these halos, and prevent its further cooling and fragmentation. 
Such close halo pairs therefore provide possible sites in which primordial gas clouds 
collapse directly into massive black holes (Mbh ~ lO'^^^M©), and subsequently grow 
into supermassive (Mbh >10^Mq) black holes by z w 6. 
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1 INTRODUCTION 

Radiative feedback plays an important role in the 
formation of stars, black holes and galaxies. Ioniza- 
tion and heating by photons with energies exceed- 
ing the hydrogen ionization threshold E > Eh af- 
fect the ability o f gas to cool and co l lapse into 
dense objects fe.g [Couchman fc ReesI Il986l: lEfsta.thioul 



19921: iThoul fc Weinberg! 1 199d : iKitavama fc Ikeuchil I2OO0I : 
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formation, the dominant coolant of primordial gas clouds 
at gas temperatures T ^10"*^ K ar e H2 r nolecules (e.g. 
ISaslaw fc Zipwl 1 19671 : iLepp fc ShuUl 1 1984 iHaiman et all 
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Il996l ). The H2 molecules can be photo-dissociated by ultra- 
violet (UV) radiation, either directly (by photons with en- 
ergies E > 14.7 eV, if the molecules are exposed to ionizing 
radiation) or as a result of electronic excitation by Lyman- 
Werner (hereafter LW) photons with energies 11.2eV <, E 
< En- In this paper we focus on this latter process, which 
operates even in gas that is self-shielded, and/or in gas that 
is shielded by a neutral intergalactic med ium (IGM) prio r 
to the completion of reionization (see Haiman et al.lll997l ), 
from radiation at i5 > Eu- 

Photodissociation feedback possibly plays an impor- 
tant role in the formation of the supermassive black holes 
(SMBHs, Mbh ~ 10^ M©) that existed at reds hift z > 6 , 
when the age of the universe was < 1 Gyr (e.g. lFanll2006l ). 
Several studies have modeled the growth of these SMBHs 
due to accretion and mergers, starting from stellar-mass 
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seed BHs left behind by the first generation of stars (e.g. 
Haiman fc Loebll200ll:lHaimanll2004i:lYoo fc Miralda-Escud^ 



2004'; Bromlcv, Somerville & Fabian '2004'; 'Shaoiro' 2005'; 
yolontcri fc Rccs 2006.; ,Li et aL. ,2007. : .Johnson fc Bromrg 
200?! ; iTanaka fc HaimaiT 20081 ). The generic conclusion is 
that it is challenging to account for the presence of 
severalx 10^ M© SMBHs by z ^6 if accretion is limited to 
the Eddington rate, unless the accretion is essentially un- 
interrupted over the Hubble time. Gravitational recoil dur- 
ing mergers can lead to the ejection of growing seed BHs 
from their par ent halo s and can exacerbate the problem 
llHaiman 2 o3; IYoo fc M iralda-Escudc 2004; Shapiro 2005; 
IVolonteri fc ReedbOod T lTanaka fc HaimanI 120081 ). An addi- 
tional difficulty is that models successfully producing the 
severalx 10^ M0 SMBHs by z >6 tend to overproduce the 
number of lower-mass BH s ([Bromley. Somerville fc FabianI 
[2004; Tanaka fc Haiman 20oi)^ 

Several alternative mechanisms have been proposed re- 
cently that evade these problems by having a period of 
rapid super-Eddington growth (e.g. Bromrti fc Loebl 20031; 



Volonteri fc Reesll2005l ;lBegelma n et al.ll2006l : ISpaans fc Silkl 
20061 ) . The reason that photodissociation feedback may pro- 



vide an interesting formation channel is that this feedback 
mechanism can affect gas coo ling and collapse in to halos 
that have Tvi^ > 10* K (e.g. lOh fc HaimanI |2002| ). With- 
out H2 molecules, gas inside these halos collapses nearly 
isothermally due to atomic cooling, with the temperature re- 
maining as high as r ~ lO'^K, which may strongly suppress 
the ability of gas to fragment into stellar mass objects (e.g. 
lOh fc Hairnan[i2002D . Instead of fragmenting, this gas c ould 
rapidly accrete onto a seed BH (|Volonteri fc Ree3l2005l ). or 
collapse dir ectly into a very massi ve (Mbh = 10'* — Iff Mp, 



black hole l|Bromm fc L ocb 200 31; iKoushiappas et al.l 



2004 



20071 ; 



Begelman et al' '2006'; Lo dato fc NataraianI |2006| . . _ _ 
Spaans fc Silk 2006; Volon teri et al.ll2008l ). "possibly with , 
intermediate state in the form of a very massive star (see 
lOmukai et al. I I2OO8I , for a more complete review). If these 
"direct-collapse black holes" indeed formed, then they would 
provide a head-start that could help explain the presence of 
SMBHs with inferred masses of several 10^ M© per comoving 
Gpc^ by 2 « 6. 

Possibly the most stringent requirement for these rapid- 
growth models is the absence of H2 molecules during cloud 
collapse. This requires the existence of a photodissociating 
bac kgrouncQ whose mean intensity exceeds J2i,lw 1? Jcrit ~ 
10^ (|Bromni fc Loebll2003l ) . where J2i,lw denotes the inten- 
sity in the LW background (in units of 10~^^ erg s~^ Hz~^ 
sr~^ cm~^, evaluated at the Lyman limit). Since LW pho- 
tons propagate nearly unobscured through the high-redshift 
(neutral) intergalactic medium until they redshift into one of 
the Lyman-series transitions of atomic hydrogen, each halo 
simultaneously 'sees' a large number of LW-sources. Each 
halo is therefore expected to be exposed to approximately 
the same LW-flux, close to the global mean level of the back- 
;round, which is lik e ly to be significant ly lower than Jcrit 
Johnson et ai]|2008l ; lOmukai et aI||2008D . 

However, some variation in the LW-fiux that irradi- 



Note that the gas could remain H2-free and close to T ~ 10*K 
if its atomic cooling rate was r educed sufficiently d ue to the trap- 
ping of Lyman line radiation jSpaans fc SilklbOOd) . 



ates individual dark matter halos is expected. For exam- 
ple, a dark matter halo with a nearby star forming galaxy 
(separated, for example, by d « lOkpc) will be exposed 
to a LW fiux exceeding the global value, and reaching 
J2i,LW ~ 10^(d/10 kpc)-^[M./(2OM0 yr"^)], where M. is 
the star formation rate in the nearby galaxy (see ii l2.1.2l be- 
low for a more detailed discussion). In the present paper, 
we model this variation and compute the probability distri- 
bution function (PDF) of the LW-fiux that irradiates high- 
redshift dark matter halos. In particular, we focus on the 
high-fiux tail of this PDF, as our goal is to answer the fol- 
lowing questions: what fraction of DM halos with Tvir > 10^ 
K, collapsing at high redshift, are exposed to a flux exceed- 
ing Jcrit? Is this fraction sufficient to account for the space 
density of severalx 10^ M0 SMBHs at 2 ~ 6 inferred from 
observations? 

In the model we present below, we account for the clus- 
tering of dark matter halos. Proper modeling of clustering 
is important when one wants to model the high LW-flux 
tail of the distribution, which is dominated by close pairs of 
halos. We also account for Poisson fluctuations in the num- 
ber of star forming galaxies surrounding a halo, and allow 
for a scatter in LW luminosity produced in a halo with a 
given mass. In the local universe, the UV luminosity can 
vary by orders of magn itude for a given stellar mass (e.g 
ISchiminovich et al.l|2007 l) , and a similar variation of the U V 
luminosity is expected for a given total halo mass. 

Spatial variations of the UV background have been in- 
vestigated previously by several authors, especially in the 
contex t of t he Ly man a forest at lower-redshift. For ex- 
ample, IZud (|l992l ) studied the expected fluctuations in the 
ionizing radiation field produced by randomly distributed 
sources, and showed that the variance is significantly in- 
creased due to radiative t ransfer effects. More recently, 
iMesinger fc Dijkstra I l|2008l ) studied the expected fluctua- 
tions of the ionizing background during the advanced stages 
of reionization using 'semi-numerical' simulations. Previ- 
ous works have studied various observable consequences of 
fluctuations in the ionizing background on th e lower red- 
shift (z <:5.5) Lyg forest (e.g. Croft 2004 ; Meiksin fc White! 
I2OO4I ; ISc"hav3l2006l ). Existing work on the evolution of the 
LW background has focused on its global (volume aver- 
aged) build-up with time, and the impact on subsequent star 
formation (e.g . Haiman et al.l I2OO0I : IMesinger et al.l I2OO6I : 
I Johnson et alBoOS ) .The main differences of our study from 



these earlier works are that (i) we include the non-linear 
clustering of sources, (ii) we specialize to compute the PDF 
of the flux as seen by high-redshift DM halos, (iii) we allow 
the UV luminosity for a given halo mass to be variable, in- 
stead of assuming a rigid one-to-one correspondence between 
halo mass and UV luminosity, (iv) we focus on the tail of the 
flux PDF, and (v) we discuss the signiflcance of this tail for 
rapid high-redshift SMBH growth (note however that the 
work presented by Mesinger fc Dijkstra (2008) did include 
effects (i) and (ii)). While this pape r was being comp leted, 
we became aware of related work bv lAhn et al.l (|2008l ). who 
studied the inhomogeneity of the LW background with cos- 
mological simulations. Their paper focuses mostly on the 
impact of these fluctuations on the star formation efficiency 
in minihalos (as opposed to the tail of the PDF and its 
significance to SMBH formation, which is the focus of our 
study). The papers seem to agree wherever there is overlap: 



the z = 10.5 LW-flux PDF shown in Fig. 11 of lAhn et al] 
l|2008f ) is consistent with our results (e.g. our Fig. [2]). 

The outhne of this paper is as follows: In § [2] we de- 
scribe our model for calculating the LW flux PDF. In § O 
we present our results. Model uncertainties and the im- 
plications of our work are discussed in § 21 before sum- 
marizing our main conclusions in § O The parameters 
for the background cosmology used throughout this pa- 
per are Oa, fit, /i, erg) = (0.27,0.73,0.042,0.70,0.82), 
consistent with 5-year data from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) (|Dunklev et al.l [20081 : 
iKomatsu et~al]l2008l ). 



2 THE MODEL 

This section, which describes our numerical modeling, is 
split into three parts. In § 12.11 we describe how we model 
the non-linear spatial clustering of sources around the po- 
sition of a dark matter halo. In § 12.21 we describe how we 
generate a discrete number of UV sources from a contin- 
uous density field, and how we assign UV-luminosities to 
individual sources by using a Monte-Carlo (MC) technique. 
In § 12.31 we discuss how we finally obtain the total intensity 
of the Lyman- Werner radiation field (Jlw) that the halo is 
exposed to. This procedure yields individual realizations for 
Jlw seen by a single halo. A distribution for Jlw is obtained 
by performing multiple Monte-Carlo realizations for a range 
of halo masses. 

2.1 Modeling the Lyman- Werner Background as 
Sampled by Dark Matter Halos 

2.1.1 The Clustering of Sources around a Halo 

Consider a halo of total (dark matter-|-gas) mass M that 
virializes at redshift z. The average number N{m,r)dmdr 
of halos within the mass range m ± dm/2 that populate a 
surrounding spherical shell of physical radius r and thickness 
dr, is given by 

A/'(m, r)dmdr = Aur^dr 

x(l + .f ^"^"^^'") dm[l + e(AJ,m,.,r)]. (1) 

Here dnsT{m, z)/dm is the Pre ss-Schechter (1974 ) mass 
function (with the modification of lSheth et al.|[2001^ ■ which 
gives the number density of halos of mass m (in units of 
comoving Mpc~'^). The factor (1 -I- z)^ converts the number 
density of halos into proper Mpc~^ . 

For simplicity, we assume a static Euclidean space in our 
calculations. Cosmological corrections affect only sources at 
distances comparable to the Hubble length from the central 
halo, and we have explicitly verified that the value of the 
global LW-background would be changed by less than a fac- 
tor of ~ 2. Since this is within the uncertainty in the flux 
attenuation due to the uncertain intergalactic H2 abundance 
(see 14. ip . we have ignored these cosmological corrections 
throughout our paper for simplicity. 

The quantity £,{M, m, z, r) denotes the two-point corre- 
lation function, which gives the excess (above random) prob- 
ability of finding a halo of mass m at a distance r from the 
central halo. In this paper we are especially interested in the 
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Figure 1. The two-point correlation function ^(M,m, z,r) is 
shown for redshift z = 10, and halo masses M = m = 1.7 X 
10* Mq. The solid (dashed) curve shows ^{M,m, z,r) obtained 
using the non-linear (linear) bias approximation, while the his- 
togram shows the two-point correlation function derived from a 
semi— numerical simulation (see text) . The non-l inear bias forma l- 
ism matches the simulation better (as shown bv llliev et al.|[2003t) . 
and therefore provides a more accurate description of clustering 
of closely separated halos. Note that in the range of r 10 — 100 
(comoving) kpc, the two-point function in the non-linear case ex- 
ceeds the linear-bias prediction by 1-2 orders of magnitude. 



high-end tail of the PDF fiux impinging upon halos of mass 
M. Since, as we demonstrate below, this tail is dominated 
by close pairs of halos, we would like to model the cluster- 
ing of closely separated halos as accurately as possible. To 
this end we use the an alytic formulation of non-linear Eule- 
rian bias developed bv llliev et aLl(j2003l l. which fits the two- 
point correlation function derived from N-body simulations 
significa ntly better than the standard linear-bias approx i- 
mation (|Mo fc White! 1 19961 : IScannapieco fc Barkanall2002l ) , 
especially at small separations. 

The benefit of using this non-linear bias formalism is 
illustrated in Figure [1] where we show, as an example, the 
two-point correlation function £,{M, m, z, r) for z = 10, and 
M = m = 1.7 X 10* M©. Here, the solid (dashed) curve 
shows the non-linear (linear) bias approximation, while the 
histogram shows the two-point corr elation function derived 
from a 'semi-numerical' simulation (iMesinger fc Furlanettd 
[2003). In this simulation, an excursion-set approach is com- 
bined with first-order Lagrangian perturbation theory to 
generate density, velocity, and halo fields at z = 10 (for 
a more detai led description of this appr oach, the reader 
is referred to IMesinger fc Furlanettoll2007^ . The simulation 
does not resolve scales below r ~ 0.14 Mpc. Clearly, the 
non-linear bias formalism matches the simulation much bet- 
ter, especially at r < 1 (comoving) Mpc. Furthermore, Fig- 
ure [1] shows that in the range of r « 10 — 100 (comoving) 
kpc, the two-point function in the non-linear case exceeds 
the lin ear-bias predict i on by 1-2 orders of magnitude. Re- 
cently, iHennawi et al.l (|2006l ) have found a large excess of 
close quasar pairs, which implies a steepening of the (pro- 
jected) two-point correlation function, which is perhaps an 
observational signature of nonlinear clustering. 

The ionizing and Lyman- Werner flux received from a 
halo of mass m at a distance r depends on the UV-emissivity 
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of this particular halo, which is likely to be variable (while 
this paper was reviewed, Lee et al. 2008 posted a paper that 
presents a detailed study of the variation of UV luminosity 
with halo mass, as we discuss next. 



2.1.2 The UV Luminosity of Dark Matter Hales 

For the UV luminosity of a dark matter halo, we adopt 
the empirical re lation -Llw ~ 8 x 10'^^ A/* erg s~^ Hz~^ 
ijKennicuttI 19981 ). where M« is the star formation rate inside 
the halo in units of Mq /yr. 

The star-formation rate in dark matter halos is highly 
uncertain (one of the reasons being radiative feedback it- 
self), especially at the redshifts of interest [z ~ 10) where 
essentially no data exists. The best constraints on the effi- 
ciency of star formation in the high-redshift Universe are 
provided by observations of z ~ 6 galaxies, which populate 
a Universe that is "only" 0.5 Gyr older than at z ~ 10. 
Among observed galaxy samples, the 2 ~ 6 population is 
therefore likely to be ar closest resembla nce to the z ~ 10 
population. Note that IStark et al.l (|2007l l have found a few 
gravitationally lensed emission-line objects, which provide 
possible candidates for 2 = 10 Lya emitting galaxies. The 
small number of candidates, and the small volume that they 
populate, prohibit strong constraints on the global star for- 
mation efficiency at z — 10. However, if the candidates are 
indeed z = 10 Lya emitters, then their abundance would 
imply that star formation is very efficient at 2 = 10 (see 
Fig 6 of Mesinger & Furlanetto, 2008), which would render 
our estimates below of the fraction of halos that "sees" an 
anomalously large Jlw conservative. 

One can reproduce the observed rest-frame UV lumi- 
nosity function of Ly man Break Galaxies at 2 ~ 6 in the 
Hubble Deep Fields (|Bouwens et all \200dt ) by simply as- 
suming that a fraction euc ~ 0.1 of all dark matter halos 
are converting /, ~ 0.1 of their gas into stars over their 
duty-cyc le of epc^hub ~ O.lth ub, where thub is the Hub- 
ble time (|Wvithe fc Loebll2006D . Furthermore, such a model 
also successfully reproduces the observed luminosity func- 
tions of Lyg emittin g galaxies at 2 = 5.7 and 2 — 6.5 (e.g 
iDiikstra et al.l |2007| ) , which suggests it provides a decent 
prescription of star formation at 2 ~ 6. If we apply the 
same star formation efficiency to the halos of interest, then 
we find that 

(W)M^2.8(^)(i±^)-^■^ (2) 

where the UV-luminosity is in units of 10^^ erg s~^ Bz'^. 

Equation ([2} implies a one-to-one relation between the 
mass m of a dark matter halo and its UV-emissivity. How- 
ever, in the local universe the UV luminosity varies by or- 
ders of magnitude for a given stellar mass (e.g. Fig. 7 of 
ISchiminovich et al.|[2007l '). 

In this paper, we model the spread in UV luminosity 
produced in halos with halos of mass m with a log-normal 
distribution, P(-Llw,26, "i). The quantity Llw, 26 here de- 
notes the LW-luminosity per unit star formation rate, in 
units of 10^^ erg s"'^ Hz~^. The probability that a halo of 
mass m has a UV luminosity in the range log Lj, ±dlog Li;/2 
('x' denotes 'LW,26') is given by 



Table 1. The relation between the parameters ctlw a^nd ^ that 
are used to model the scatter in UV luminosity for a given halo 
mass (see eg. [3ll. 



0"LW 


0.0 


0.25 


0.50 


0.66 




1.0 


0.88 


0.38 


0.21 



P(logLa;, m)dlogLa; = 

rflog |- -(logL^ -log[.F X (L^)])'' ] 

where log denotes logio- In this paper, our fiducial model 
adopts (Jlw = 0.50. Hence, the Lyman-Werner fiux is 
larger/smaller by an order of magnitude than the value given 
by equation ([2| for 2.5% of the halos (for all m). The fac- 
tor reduces the mean luminosity (L^c) relative to the no- 
scatter prediction in equation (2] Since the inclusion of scat- 
ter increases the predicted abundance at fixed UV luminos- 
ity, this factor is necessary to maintain consistency between 
our model a nd the observed 2 = 6 rest-frame UV luminos- 
ity function l|Bouwens et al.ir2006l ). In our fiducial model, we 
find J- = 0.38; increasing ctlw requires lower values of .7^ (see 
Table [l|. We emphasize that the distribution given by equa- 
tion[2]is clearly ad-hoc - it is meant to crudely represent the 
luminosity scatter one may reasonably expect among high- 
redshift halos. The impact of varying the parameter ctlw is 
discussed in § 13.31 

2.2 Monte-Carlo Realizations of the Model 

The environment of a halo is sampled by A'^,- concentric 
spherical shells from r = rmin out to a maximum radius 
''max- We denote the radius and thickness of shell number 
j by rj and drj, respectively. Furthermore, the mass func- 
tion dnsT{m, z) / dm is sampled by Nm mass bins that are 
spaced evenly in logm. Mass bin number i contains halos in 
the mass range logm^ ± dlogmi/2. 

2.2.1 The Clustering of Sources around a Halo 

Equation Q gives the average number of halos 
N{m,r)dmdr within the mass range m ± dm/2 and 
within the spherical shell at r ± dr /2. The actual number 
N{mi,rj) in shell number j in the ith mass bin fluctuates 
around this average, and in the Monte-Carlo simulation we 
populate this shell with halos by drawing from a Poisson 
distribution with a mean {N)i,j = N{mi,rj)dmidrj. 

2.2.2 The UV Luminosity of the Surrounding Halos 

As mentioned above, only a fraction edc ~ 0.1 of all halos 
are actively forming stars at a given time. We determine 
whether a halo is "on" or "off" by generating a random 
number Q ^ R ^ 1. \i R ^ edc, then the halo is "on" (i.e. 
star forming) and vice versa. For each star-forming halo 
that occupies radial shell number j and mass bin number 
i, we generate UV luminosities Lxinrii) from the probabil- 
ity distribution given in equation ((3} above. Specifically, we 
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convert uniformly distributed random numbers ^ i? ^ 1 
to luminosities log La; (mi) using 



Table 2. Fiducial model parameters (see ^ l3.1l for discussion). 



R = 



log Lx {mi) 



P{u, 7ni)du. 



(4) 



This equation can be inverted analytically and yields 
\ogLx{mi) = log[:'^x (L^)(mi)] + cr^\/2erf"^(2i?-l), where 
erf~^ is the inverse error function, and {Lx}{m) is given by 
equation ([2}. 



2.3 The UV-flux at the Halo 



Using the prescription of § 12.21 for populating the environ- 
ment of a (single) halo with UV-emitting galaxies, we can 
calculate the Lyman- Werner flux that is "seen" by the halo 
from 



Jlw — -. — 



EE 



JV(ti 



E 



X A{R- 



eDC,k) 



(5) 



Here A(a;) is the Heavyside step function, and the factor 
l/47r preceding the summation symbols converts the units 
of Jlw into erg s~^ cm~^ Hz~^ sr~^ (the units of Llw 
are erg s~^ Hz~^). The sum over k reflects the fact that 
any given mass-radius bin may contain > 1 halos. For suffi- 
ciently large radii and low masses, we often have {N)i^j ^ 1. 
In this case, we do not generate UV luminosities for indi- 
vidual halos. Instead, whenever {N)i^j ^ 10, the total UV 
luminosity from that shell is given by its mean value, i.e. 
Lx,i,j = {N)i^j X EDO X {L^){mi). This approximation speeds 
up the calculation, and does not affect the Jlw PDF in the 
tail. The main reason is simple: the high- Jlw tail of the 
PDF arises from Poisson fluctuations in the number of ha- 
los in the inner radial shells, and from fluctuations in the 
UV-luminosity of these sources (see below). 

Equation (|5]) gives the Lyman- Werner flux that is seen 
by a single halo. We repeat the Monte-Carlo calculation A^mc 
times in order to derive an accurate PDF of Jlw- We would 
like to investigate whether halos exist that are irradiated 
by a LW-background that exceeds Jcrit ~ W^, so that the 
gas inside these halos may collapse directly into a massive 
black hole. This scenario was introduced as a possible for- 
mation mechanism for the seeds of z — 6 supermassive black 
holes (see § [!}. These black holes are extremely rare, with 
a number density of ~ 1 (comoving) Gpc~^ (e.g. Fan et al. 
2001, 2003). Hence (approximately) only one halo per cGpc'^ 
needs to yield a seed. For comparison, there are ~ 10^ ha- 
los at 2 = 10 with virial temperatures exceeding Tvir = 10'' 
K in a 1 cGpc^ volume. Ideally, we would therefore want 
to perform A'mc ~ 10^ Monte-Carlo simulations, and ver- 
ify whether at least one halo is exposed to a LW-flux that 
exceeds Jcrit. However, this is computationally prohibitively 
expensive. In practice, we instead utilize A^'mc ~ 10^ real- 
izations, and complement our Monte-Carlo simulations with 
analytic calculations as described in § 13.31 



3 RESULTS 

We use Nm = 500 and A^'^ = 100 throughout this paper. We 
have verified that our results are not sensitive to the precise 
choices of Nm and A^^. 
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Figure 2. The probability distribution (PDF) of the Lyman- 
Werner flux Jlw (in units of J21 lw = 10~^^ erg cm~^ 
sr~^ Hz~^) as seen by halos of mass M = 4 X lO'^ Mq at z = 
10. The black solid (red dotted) histogram corresponds to our 
fiducial model with (without) clustering. The dashed vertical line 
at J2i,LW = 40 denotes the mean value of the LW background 
(see text). The figure shows that (i) the vast majority of halos 
see a LW fiux that is within a factor of 2 of the mean value, (ii) 
the clustering of halos boosts the tail of the flux PDF by more 
than an order of magnitude, and {iii) an exceedingly small, but 
non-zero fraction (~ few xlO"'^) of dark matter halos see a LW 
flux boosted to levels exceeding J21 lw > Jcrit = 10^. 



3.1 The Model Parameters 

In our fiducial model, we investigate the flux PDF seen by 
halos of mass M = Mh = 4 X 10^ M© at 2 = 10. These halos 
are just massive enough to excite atomic H cooling processes 
(i.e. Tvir = 10*K, e.g. eq. 25 in Barkana & Loeb 2001, for a 
mean molecular weight /i = 1.2). As mentioned above, the 
number density of sources more massive than Mh at 2: = 10 
is ~ 1 cMpc""^. The impact of varying M is investigated in 
§E3] 

The minimum distance of the center of the halo to the 
center of another nearby halo, rmin, has to be specifled in our 
model. In the absence of peculiar velocities, one may expect 
''min > '"vir.M. Howevcr, closc palrs of halos can merge, and 
share a common halo, making rmin < ^vir.M possible. On 
the other hand, (very) close halo pairs are short-lived, and 
will merge on a time scale shorter than eociiiub. The pre- 
cise choice of rmin is therefore somewhat arbitrary; in the 
flducial model we adopt rmin = 2rvir ~ 2.0 kpc, and we will 
investigate the impact of varying rmin in § 13.31 

In our calculation, we only consider halos out to rmax = 
(Alw,i — ^0)c/[^pH{z)] ~ 18 pMpc, in which Alw,i = 
1110 A is the wavelength corresponding to the lowest en- 
ergy Lyman- Werner transition. Lyman- Werner photons that 
were emitted at r > rmax have redshifted into one of the 
Lyman-series resonances prior to reaching the halo, and are 
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Figure 3. Two realizations of the environment of a fiducial Af = 4 X 10^ Mq fialo at 2 = 10. The left {right) panel shows a halo that 
"sees" a J21 lw = 10^ (J21 lw = 40). The fiducial halo is represented by the filled square sA, x = y = kpc. Open circles denote 
surrounding galaxies that are not actively forming stars, and therefore not contributing to the local LW flux. Solid circles denote halos 
that are actively forming stars; their UV luminosity is represented by the size of the circle (the largest circles are the most luminous, as 
illustrated by the graphic scale above the panels). The halo that sees J2i,lw ~ 10"^ has a nearby (r ~ 3 kpc) neighbor that is luminous 
(i'LW,26 = 10^)1 while the halo in the right panel has no luminous nearby neighbor, and sees J21,lw ~ 40. The latter configuration is 
most common. 



converted efficiently into photons with energies < 10.2 eV 
before reaching our halo. In reality, 18 pMpc is an upper 
limit, and rmax sho uld be a strong function of frequency 
l|Haiman et al.|[l997l ). The impact of LW radiative transfer 
on our results is discussed further in § 14.11 in practice, the 
high-flux tail of the PDF is insensitive to the choice of rmax. 

Finally, we assume that the minimum mass for star 
forming halos is mmin = 4 x 10^ M©. This model there- 
fore represents a universe in which no star-formation occurs 
in "minihalos" with Tvir = 10^ K (e.g. due to the presence of 
the LW background itself or due to an X-ray background). 
The impact of varying rrimin is investigated in § 13.31 The pa- 
rameters of our fiducial model are summarized in Tabled 

3.2 The Flux PDF in the Fiducial Model 

In Figure [H we show our main results: histograms showing 
the fraction of Monte-Carlo reafizations that yield a Lyman- 
Werner fiux Jlw as a function of Jlw- The solid histogram 
shows the fiducial model. The distribution in Jlw peaks at 
around J2i,lw = 40, and 95% of the models lie in the range 
35 jjJ2i,LW ^45. The distribution is highly asymmetric with 
a long tail towards high values of J2i,lw. 

The sharp cut-off on the 1ow-J2i,lw side is due to the 
fact that the IGM is assumed to be transparent to LW- 
photons, which causes a large number of distant sources to 
contribute to the intensity of the local LW radiation field 
(for reference, the total number of halos with m > mmin 
within a sphere of radius rmax = 18 pMpc is ^10^!). Poisson 
fluctuations in a large number of distant sources, as well as 
variations in their UV luminosity, only result in a small vari- 
ation in Jlw. Indeed, the dashed vertical line at J2i,lw = 40 
corresponds to the value of Jlw in the absence of any scatter 
in UV luminosity of individual halos (i.e. (tlw = 0.0), and 
with no Poisson fluctuations in the number of halos. This 
can be thought of as the level of the global mean LW back- 
ground. No halos see a flux that is significantly smaller than 
this mean background value of J2i,lw = 40 (i.e. among our 



10^ realizations, and with our flux-bin resolution, there are 
no instances when the flux falls below the lowest flux-bin. 
Note that this value differs for the runs that do and do not 
include clustering). 

On the other hand, the long tail towards high-J2i,LW 
arises from Poisson fluctuations in the number of close ha- 
los, and from the assumed scatter in their UV luminosity. 
The highest values of J2i,lw are seen by halos that have a 
luminous nearby neighbor. This is illustrated graphically in 
the left panel of Figure [3] where we show the environment 
of a halo that is irradiated by a LW-flux close to the maxi- 
mum value we find, J2i,LW ~ 10^. The smaU (red) square at 
X — y = kpc denotes the position of the central halo, while 
open circles denote halos in which no star formation occurs, 
and are therefore not producing any LW fiux. The filled cir- 
cles denote halos that are forming stars, with the magnitude 
of the UV-luminosity represented by the size of the circles. 
For reference, the sizes corresponding to I/lw,26 = 1.0, 10.0 
and 100.0 are shown as labeled above the figure panels. 
The central halo is illuminated by a bright (Llw,26 ~ lO'^), 
nearby (r « 3 kpc) star forming halo0 For comparison, the 
right panel shows the environment of a much more common 
halo that sees Jlw, 21 = 40. Here, the nearest source is lo- 
cated at a distance of r = 25 kpc and has a UV luminosity 
of J/lw,26 = 10 (this source only contributes ~ 2% of the 
total LW fiux seen by the central halo at r = 0). 

Note that this figure also illustrates the importance of 
clustering. Figure |3] shows the distribution of halos in a sur- 
rounding volume that is ~ 0.35 comoving Mpc^; without 
clustering, according to the Sheth-Tormen mass function, 

^ Our spherically symmetric clustering model generates only ra- 
dial distances; we have artificially generated random azimuthal 
angles to place the halos in the x — y plane in Figure |3] Hence, 
Figure |3] does not represent an accurate projection or slice of a 
realistic 3D source distribution; it is meant only to graphically 
illustrate the difference in the environment of halos illuminated 
by a high versus a low Jlw- 
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this volume would contain, on average, only ~ 1.4 dark mat- 
ter halos with Tvir ^lO^K. In the presence of clustering, we 
find several tens of halos in both realizations shown in the 
two panels. 

Figure [2] shows that even though the distribution of 
Jlw drops off rapidly with increasing Jlw, the fraction few 
X 10~^ of all halos see J2i,lw ^ 10'^. This corresponds to 10^ 
halos in a volume of 1 Gpc"^. Note that because the high- 
Jlw tail is dominated by close pairs of halos, accounting for 
clustering is quite important in the present calculation. This 
is quantified by the red dotted curve which shows the fiducial 
model, but without source clustering (i.e. ^{M,m, z,r) — 
in eq. In the absence of clustering the fraction of halos 
that sees a given Jlw above the mean value is typically 
smaller by about 1.5 orders of magnitud^f]. 

Our main result, presented in this section, is that a 
non-negligible fraction of halos, i.e. a few xlO~^, may see a 
local LW fiux that is higher by almost two orders of magni- 
tude compared to the mean global background (correspond- 
ing to a model in which Poisson fluctuations in the num- 
ber of sources, their clustering, and variations in their UV- 
luminosity are ignored). In the next section, we investigate 
the sensitivity of this result to the assumed model parame- 
ters. 

3.3 The Impact of Model-Parameter Variations 

The impact of varying the model parameters on the very 
high Jlw tail of the fiux distribution can be efficiently in- 
vestigated using an analytic approach. This is because the 
high- Jlw tail is dominated by halos with a nearby star- 
forming halo, and, in this limit, the Poisson probability of 
having two or more such nearby halos is negligibly small. 
The probability that a halo has at least one neighbor that 
would alone expose it to a Lyman- Werner fiux in the range 
log Jlw ± dlog Jlw/2, is given by 

J'i(JLw)dlog Jlw = dlogJLwxeocX 

JV„, JV, 

E E (iV)i,je-<^>-P(logL[JLw,r,],L.K]) (6) 

where (A'')i,j — J\f{mi,rj)dmidrj (see I2.2.ip . r(jmin) = 
' min , m(jlmin) and P{\ogL[Juiv,r],Lx[m]) dlog 

Jlw denotes the probability that a halo of mass m that 
is located at a distance r produces a LW luminosity in the 
range log Jlw ± ci log Jlw/2. This halo must have a Lyman- 
Werner luminosity that is L = IGvr^r^ Jlw (in erg s~^ Hz^^), 
and P{logL[Ji,w,r],Lj:[m\) dlog Jlw = P(logL, m)dlogL, 
for which the expression was given in equation (|3]). 

In Figure |31 we compare the probability 
A (Jlw) dlog Jlw obtained from equation © {dotted 

^ The enhancement of the high-flux tail of the PDF due to 
clustering varies only weakly with J21 lw- The reason for this 
weak dependence is that especially the tail-end of the PDF at 
Jlw, 21 ^j200 arises mainly from single UV-sources that are in 
close proximity (2 kpc< r < 7 kpc, see ii l3.3l l to the central halo. 
Over this narrow range of separations, the two-point function 
^(M,m, z,r) — and hence the expected number of surrounding 
halos - varies relatively weakly (the clustering length for halos of 
mass M is ~ 300 pkpc, see Fig. [TJ 




100 1000 

Figure 4. The probability Pi (Jlw) for a halo to have a neigh- 
bor that would alone expose it to a Lyman- Werner flux of Jlw is 
shown by the dotted curves (obtained from eq. |6j , and compared 
to the flux PDF obtained from Monte-Carlo simulations (his- 
tograms). The flducial model is shown in the middle panel, while 
in the top panel {bottom panel) the luminosity scatter ctlw is in- 
creased (decreased) to ctlw = 0.66 (ctlw = 0.25). Both methods 
yield consistent results for Jlw, 21 ^500 for all models, demon- 
strating that the high- J2i_lw tail of the Jlw distribution is dom- 
inated by a single nearby star forming halo, and that the analytic 
approach can be used to extrapolate the probability distribution 
of Jlw beyond J21,lw ^500- The dashed curve in the middle 
panel shows a variation of the analytical prediction that accounts 
for pairs of neighboring halos producing the combined flux Jlw- 
The agreement at lower values of J21.LW is better, indicating 
that at Jlw, 21 ^500, often two halos contribute signiflcantly to 
the total flux. 



curves) with that obtained from the Monte-Carlo simula- 
tions {histograms) for the fiducial model {middle panel), 
and for a model in which ctlw is changed to ctlw ~ 0.66 
(ctlw = 0.25) in the top {bottom) panel. Figure |4] shows that 
the two methods yield consistent results for Jlw, 21 ^500 
in all three cases. This agreement demonstrates that the 
high fiuxes must be dominated by a single nearby source, 
and that we can safely use equation ([6]) to compute the 
JLW-probability distribution for J2i,lw <;500. At lower 
values of J2i,lw, the Poisson probability of having more 
than one nearby UV-bright star forming halo can become 
non-negligible. To examine this possibility, we define the 
probability that a halo has either at least one neighbor 
that would alone expose it to a Lyman- Werner fiux in 
the range log Jlw ± dlog Jlw/2 (as before), or it has at 
least one pair of neighbors whose combined fiux equals 
log Jlw ± dlog Jlw/2. In the latter case, the two halos are 
required to produce intensities of yJhw and [1 — y] Jlw, re- 
spectively, with < y < 1. For a given fractional allocation 
y of the fiux between the two halos, P2( Jlw) dlog Jlw is 
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Figure 5. The cumulative probabilities for the mass {top panel) 
and location (bottom panel) of the single halo that produces 
•^21 LW = 10^ in the fiducial model (solid curves) and for a model 
in which rmin = 20 kpc (dashed curves ). The light grey regions 
denote the central 68%, or Itr, confidence interval. The figure 
shows that in the fiducial model, the most probable cases are 
low-mass halos (in the range 3 X 10® < m < 6 X IO^Mq) lo- 
cated very nearby (2 < r/kpc < 7). The halo mass and distance 
are both boosted for the model in which r^[^ is increased to 20 
kpc. 



given by P2(i/) oc {Pi(j/Jlw) x Pi([1 - 2/]Jlw)}, and the 
total probability for at least one such pair of neighbors is 
given by P2,tot = Jq dyP2{y). As an example, the dashed 
curve in the middle panel of Figure |4] shows the probability 
Pi + P2,tot as a function of Jlw- Clearly, accounting for 
pairs of halos provides a better fit to the MC-simulations 
at lower values of J2i,lw. 

In Figure O we investigate the most likely host halo 
mass and distance of the nearby UV source producing the 
critical flux. In the top panel, we plot the cumulative proba- 
bilities for the mass of the halo that produces J2i,lw = 10^ 
in the fiducial model (solid curve), and for a model in which 
Tmin = 20 kpc (dashed curve). This cumulative probabil- 
ity is given by replacing m,nin — > m as the lower limit in 
the mass integral in equation ©. The solid curve shows 
that in the fiducial model, the halo mass lies in the range 
3 X 10^ Mq < m < 6 X 10^ Mq, while its most likely position 
is 2 < r/kpc < 7 (the shaded ranges enclose 68% of the 
total probability around p(> x) = 0.50). The most likely 
halo mass is increased to lie in the range ~ 10^° Mq < m < 
8 X 10^*^ M0 for the model in which rmin = 20 kpc. 

In the bottom panel, we plot the cumulative distribution 
that the star forming halo was located at a distance that 
exceeds r (given by replacing rmin ^ r as the lower limit 
in the radial integral in eq. |6}. In the fiducial model, the 
distance to the dominant star forming halo lies between 2 < 
r/kpc < 7, while this range is increased to 20 < r/kpc < 
40 < for the model in which rmin ~ 20 kpc (68%-ranges are 
shaded, as above). 
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Figure 6. The impact of varying the model parameters, M (top 
panel), cTlw (second panel from the top), rmin (third panel), and 
"^min (fourth panel) on the Jlw probability distribution is shown. 
This figure demonstrates that our results are most sensitive to the 
precise values of Ad and rmin • 

Increasing M (rmin) significantly 
boosts (reduces) the high-flux tail. However, even an increase of 
rmin by a factor of 10 results in a non-negligible fraction of halos 
that are exposed to a J21 lw = i^crit = 10^ (see text). 



It is interesting to point out that in the fiducial model, 
the parameter combination r = 3 kpc and m = lO^M© (for 
which p(> r) = p(> m) = 0.5), results in J2i,lw ~ 200. 
Hence, in order for such a halo to produce a local flux of 
J2i,hw ~ 10^, it must be ~ 5 times brighter than average. 
Since this corresponds to ~ 1.4(jlw, this implies we are not 
very sensitive to the tail-end of the assumed log-normal dis- 
tribution of the UV-luminosity of star forming halos (indeed, 
in Fig. [S] we explicitly show that the tail of the flux PDF is 
not sensitive to (Jlw). 

In Figure |6l we explore the impact of varying the model 
parameters, M (top panel), ctlw (second panel from the top), 
as shown above in Fig. rmin (third panel), and rrimin 
(fourth panel). The parameter dependencies seen in this fig- 
ure are easily understood. 

Figure |S] shows that: (i) the tail end of the flux PDF is 
boosted slightly for M = 8 x lO^M© (dotted curve) and sig- 
nificantly for M = 10^ M© (dashed curve), which reflects the 
fact that sources are clustered more strongly around more 
massive halos; (ii) increasing ctlw enhances the tail of the 
flux PDF, but only slightly. The enhancement arises because 
increasing ctlw enhances the probability of having a nearb' 



roy 



halo that is (significantly) brighter than average in the U 



^ Consider a halo of mass m = 4 X W"^ Mq that is located at 
r = 2 kpc. In the fiducial model, this source needs to be ^60 times 
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Figure 7. We show the LW flux PDF for the fiducial model 
{solid black histogram, as in Fig. [2]l, and for a model in which 
intergalactic H2 reduces the photodissociation rate due to far- 
away sources, which in turn reduces the value of the global mean 
LW background by about an order of magnitude. However, the 
high- tail of the flux PDF depends entirely on the statistics 
and nature of very nearby star-forming halos, and LW— radiative 
transfer effects from the IGM do not affect this tail of the Jlw 
probability distribution at all. 

(tii) increasing rmin decreases the tail at high Jlw- This is 
because increasing rmin eliminates the contribution of very 
nearby halos to the Jlw distribution (this conclusion can 
already be inferred from Fig. [5] which shows that the most 
nearby halos dominate the single-halo flux); (iv) increasing 
reduces the fiux PDF at J2i,lw ;$500, but the high- 
flux tail is not significantly affected. This is because the 
highest values of J2i,lw require Lyman- Werner luminosi- 
ties that are L = 6 X 10^'^(rmin/2 kpc)^(J2i,Lw/10^) erg s"^ 
Hz~^. Hence, producing these fluxes requires the presence of 
a source that is forming stars at a minimum rate of ~ 0.7Mq 
yr~^. Prolonged star formation rates of this magnitude can 
only occur in halos that are more massive than 2 x 10* M©. 
Therefore, any variation of mmin has a negligible impact on 
the flux PDF at J2i,lw ^500, as long as mmin 
(which seems likely prior to the completion of reionization, 
e.g. Mesinger & Dijkstra 2008). 



4 DISCUSSION 

4.1 Lyman- Werner Radiative Transfer 

The dissociation of H2 molecules follows photoexcitation 
of their electronic states. There are 76 of these Lyman- 
Werner transitions, the majority of them having energies 
that lie blueward of the Ly/3 resonance (see, e.g.. Fig. 1 of 

brighter than average to produce Jlw, 21 <;10^i which corresponds 
to ^3.6(T. On the other hand, this same source needs to be ~ 100 
times brighter than average in the model in which (Jlw = 0.66 
(because the factor J^ scales the mean UV luminosity to remain 
consistent with the observed z = 6 luminosity function of LBGs, 
see Table[T]l, which in this model corresponds to ^3. la. Boosting 
f LW therefore boosts the number of sources capable of producing 
Jlw,21 ;^10^ (but only slightly). 



iHaiman et"aLll200Cl ). The majority of photons that are ca- 
pable of dissociating H2 molecules have therefore redshifted 
into the Ly/3 resonance over a distance that is much smaller 
than the rmax that was used in our model. In fact, rmax 
should have a strong frequency dependence (for example, 
some of the LW lines are just on the red side of an atomic 
H Lyman line; photons could not have suffered significant 
redshift before reaching such frequencies without scattering, 
implying small rmax at the frequency of these LW lines). 

The overall conclusion is that the total photodissocia- 
tion rate due to far-away sources is reduced, and that our 
simple modeling which ignores radiate transfer in the IGM 
overestimates the value of the global mean LW-background. 
Additionally, a trace fraction of molecular hydrogen in the 
IGM (a;H2 ~ 2 x 10"^) is sufficient to produc e an effective 
optic al depth in the LW lines of ~ 1 — 2 (jRicotti et al.l 
I2OOII ). Furthermore, X112 may be enhanced inside "relic H2 
regions", which increases the volume-a veraged effective op- 
tical depth of the IGM in the LW lines (I Johnson et al.ll2007 . 
also see Shapiro & Kang 1987). iKuhlen fc Madaul poosf ) 
have shown recently that X-rays produced by the first gen- 
eration of accreting black holes may boost the fraction of 
intergalactic H2 by a factor of ~ 20 around such "mini- 
quasars". If the IGM is optically thick in the LW lines, then 
rmax becomes an even stronger function of frequency. 

We next investigate the assumption that the IGM is 
optically thick in every one of the H2 lines. Note that 
this is a very conservative assumption, requiring that the 
photo-dissociated intergalactic H2 is continuously replen- 
ished, which is unlikely to be the case (Haiman et al. 2000); 
we purse this assumption here because it maximizes the im- 
pact of RT effects. 

In Figure[7]the red-dotted histogram shows the LW-flux 
PDF that is obtained when taking into account this strong 
frequency dependence of rmax- We modeled this frequency 
dependence by computing the maximum distance rmax(fc) 
from which LW[k] photons (the k*'' LW transition, when the 
76 LW lines are ranked simply in order of increasing energy) 
can be received, that is, Hrmax{k) / c — (Afc — Afc+i)/Afc-|-i. 
Here, Xk denotes the wavelength of the k**^ LW transition. 
LW[k] photons cannot by received from sources at larger 
radii, as these photons redshifted through the LW[k-|-l] 
resonance first, which makes the IGM effectively optically 
thick to sources beyond rmax(fc). A source at radius r then 
contributes the following flux to the local LW-background: 
Jlw = 4^ ^^^2 fc*^'"" 1 where fctot is the total photodissocia- 
tion rate as summed over all 76-transitions, and where fcdiss 
is the photodissociation rate that one obtains if one sums 
over only over those lines for which r < rmax(fc)- 

Figure [7] shows that in the above, conservative, opti- 
cally thick limit, our neglect of LW radiative transfer over- 
estimates the value of the mean LW-background by about 
an order of magnitude, and underestimates the width of the 
flux PDF. However, the high- Jlw tail of the flux PDF de- 
pends entirely on the statistics and nature of very nearby 
star-forming halos (within a few to a few tens of kpc) . LW- 
radiative transfer effects from the IGM do not affect this 
tail of the Jlw probability distribution at all, even in this 
extreme limit. 

We also note that LW photons dominate the H2- 
dissociation rate when the medium is optically thick and 
lacks radiation above 13.6eV. Indeed, the critical value of 
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Jcvit ~ 10'^ corresponds to the photo-dissociation time due 
to LW radiation being shorter than the free-fall time of 
the gas. In our case, the column density of neutral HI 
out to the nearby star-forming neighbor is likely negli- 
gible, and direct photo-dissociation by photons with en- 
ergies E > 15.4 eV will boost the total Il2-dissociation 
rate (the latter dominates for neutr al HI column densities 
A''hi <10"cm~^: lHaiman et al.lll997l ). This may help to keep 
the halo gas _ff2-free, at least during the early stages of its 
collapse (until contraction results in a column density ex- 
ceeding iVm ~ lO^^cm"^). 

Lastly, it is worth pointing out that since the tail of 
the flux PDF is due to a single neighbor, the impinging 
LW radiation field is highly anisotropic. The impact of an 
anisotropic radiation field on the gas could be quite different 
from that of an isotropic background (e.g., perhaps the 'far 
side' can be shielded more efficiently against the impinging 
photodissociating flux, which could allow gas to cool and 
fragment) , but a quantitative assessment of this 3D behavior 
is beyond the scope of this paper. 

4.2 The Evolutionary Phase of the Halo/Timing 
Issues 

We have found that a fraction of a few x 10^^ of all z — IQ 
dark matter halos with Tvir ^ 10^ K are irradiated by a LW 
flux that exceeds Jcrit (Fig.|S]). However, the direct-collapse 
black hole scenario requires that the gas is irradiated and 
kept H2-free for the entire ~free-fall time while it collapses 
to the center of the dark matter halo. This additional re- 
quirement reduces the fraction of dark matter halos in which 
these black holes may form. 

To quantify this reduction, we identify three relevant 
time-scales: (i) the duration over which halos are produc- 
ing UV luminosity (tium). In our model, this corresponds 
to tium = EDcihub; (ii) the duration over which the irradi- 
ation is needed to avoid fragmentation (tfrag); and (iii) the 
spread in formation times of neighboring halos (tform). In 
our model, we implicitly assumed that tform = ihub (since 
we only assumed that a fraction edc of all halos are actively 
forming stars). The overall likelihood for a combination of 
two halos in which one prohibits fragmentation in the other 
is given by the fraction of occasions in which one halo formed 
within a time ti^-m — ifrag prior to the other during a time 
iform. If the formation times of both halos are uncorrelated, 
then this probability is given by g = (tium — tfrag)/tform. If 
we write ffrag = iium//, then we find that q = (1 — 1//). 
In this case, the reduction due to the synchronization re- 
quirement is small for / ^2. Furthermore, the reduction is 
significantly less, if the formation times of the two halos are 
correlated (which is very likely, since both the collapsing 
and the star-forming halo formed in the same large-scale 
overdense region of the Universe). 

4.3 The Importance of Photoheating 

The halos of interest that are irradiated by a large LW-fiux 
are also illuminated by ionizing radiation, with an ionizing 
flux that is lower by a factor of /esc (the escape fraction of 
ionizing photons) t han Jlw- Th i s esca pe fraction is highly 
uncertain: recently, IChen et al.l l|2007h have derived a dis- 
tribution function for /esc from the cumulative Hl-column 



density distribution along sightlines to (long-duration) 7-ray 
burst (GRB) afterglows. Long-duration GRBs are thought 
to be associated with massive (M ^10 — 3OM0) young 
stars, and are therefore thought to occur within star form- 
ing regions of galaxies (see e.g. Woosley & Bloom 2006). 
GRB afterglows may therefore probe the HI column-density 
along sightl i nes to the production sites of ionizing photon^- 
IChen et al.l |20o3) show that /esc may obtain a wide range of 
values, and that (/esc) = 0.02, in good agreement with theo- 
retical calculations of the escape fraction of ionizing photons 
from galaxies in adaptive mesh refin ement hydrodynamical 
simulation^ by iGnedin et al.l (|2008D . 

From the above we conclude that the halos of interest 
may be irradiated by a ionizing flux that likely has a broad 
distribution between Jion,2i = — 10'', with (Jion,2i) ~ 20. 
To investigate the impact of this photoionizing flux on the 
dynamics of the gas in collapsing halo, we performed a 1- 
D hydrodynamica l simulation with a code that was orig- 
inally written by iThoul fc Weinberg! (|l995t ) and modified 
as described in Dijkstra et al. (2004, where it was used 
to study the impact of photoionization feedback on the 
formation of dwarf galaxies). We illuminated a gas cloud 
of mass M = 10* with an ionizing flux having am- 
plitudes of Jion,2i = 10 and Jion,2i = 10^, and spectral 
index a = 1 (i.e. the ionizing flux density is given by 
J{v) = J2i,ion(!^/«^L)"" X 10-2' erg S-' Hz-i cm-2 sr"'). 
This ionizing flux was switched on at a time = eocihub 
before the cloud would have collapsed to r = in the ab- 
sence of pressure. 

The result of these calculations are shown in Figure [S] 
Solid curves show the time evolution of the radii of gas shells 
enclosing 30%, 40%, 50% and 100% of the total halo mass 
in the absence of an externally generated ionizing radiation 
field. The dashed {dotted curves) represent the models in 
which Jion,21 = 10 (Jion,2i = 10''). The fi gure shows that 
even when Jion,2i = 10 (Jion,2i = 10^), stiU 50% (40%) of 
the gas was able to collapse, which implies that photoion- 
ization feedback does not prevent the gas from collapsing, 
even with extremely high UV fiuxes. The main reason for 
this is that photoionization feedback is most efficient when 
it can operate for a prolonged time: e.g. in our spherical col- 
lapse model, this feedback mechanism has the largest effect 
when the radiation field is switched 'on' when (or before) the 
gas completely decouples from the Hubble fiow and starts 
to collapse. However, in the final eocihub of the cloud's col- 
lapse, the gas in the halo has already collapsed to densities 
that greatly exceed the mean cosmic gas density in the Uni- 
verse, and cooling processes are more efficient. Furthermore, 
at these high gas densities self-shielding, which is ignored in 
our simulation, will become important which would make 
photoionization feedback even less important. 

^ Furthermore, evidence exists that the majority of GRBs occur 
in sub-L* galaxies (see Chen et al. 2007, and references therein), 
which are likely to bear a closer resemblance to the z = 10 star- 
forming galaxies in our model, than the massive LBGs for which 
escape fractions of ionizing photons have a lso been determined on 
t he basis of thei r obser ved Lyman Break l lShaplev et al.ll2006l) . 
^ IWhalen et"all l l2004h found /esc ~ 1 for galaxies at 2 ~ 20, 
using one-dimensional hydrodynamic calculation. However, their 
calculations focused on the first UV-sources in the Universe which 
populated much less massive "minihalos". 
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Figure 8. Solid curves show the time evolution of the radii of 
gas shells during the collapse of 10* Mq halo, in the absence of 
an externally generated ionizing radiation field (i.e Jion,2l = 0). 
Each set of curves corresponds to shells (initially) enclosing 30%, 
40%, 50% and 100% of the total halo mass. All gas in the halo has 
collapsed hj z = 10, when tjmb ~ 500 Myr. The dashed {dotted 
curves) represent models in which an ionizing background, whose 
amplitude was Jion,2i = 10 (Jion,2l = lO'^), was switched at 
a time At = eDC^hub before the cloud would have collapsed to 
r = in the absence of pressure. The figure shows that even when 
^ion,2i = 10 {Jion,2i = 10^), Still 50% (40%) of the gas was able 
to collapse. 



Lastly, it is wortii pointing out tiiat lialos form, and 
cluster strongly, at intersections of filaments. The dominat- 
ing LW-flux will therefore be emitted by halos along the 
filaments. On the other hand, their ionizing flux will most 
likely escape into the low density regions, perpendicular 
to the filaments. This argues for a value of Jion,2i/J2i,LW 
that is lower than (/esc). Furthermore, the simulations by 
iGnedin et al.l (|2008l ) suggest that the escape fraction of ion- 
izing photons varies significantly from sightline-to-sightline, 
with a small fraction of sightlines having fasc ~ 1, and a 
significant fraction having /osc ^10^'^ (mainly because of 
intervening HI and Hel). This also implies that it is likely 

that Jion,2l/ J21,LW ^(/esc). 



4.4 The Transition from Massive to Supermassive 
Black Hole 

Even if a lO" - 10® Mq black hole forms as a direct re- 
sult of the collapse of a gas cloud that was irradiated by 
a LW flux that exceeded Jcrit, the black hole still needs to 
grow to 10® M0. This requires the black hole to consume a 
tremendous amount of gas over ~ 0.5 Gyr (the total time 
that elapsed between z = 10 and z = 6). This implies that 
the original host halo in which the black hole formed, must 
merge with other halos into a massive 10^^ - 10^^ Mq host 
halo (which also agrees with the mass of the host halos of the 
2 ~ 6 quasars; see, e.g., the review by Haiman & Quataert 
2004). Interestingly, we found that those M = 10* halos 
that are exposed to a LW fiux exceeding Jcrit, are the ones 
with a nearby star forming galaxy. The halo mass of this 
nearby galaxy is possibly as large as 10^^ Mq (see § 13. 3|) . 



and the 10 — 10 M© black hole is likely to end up inside 

this larger da rk matter halo . 

Utilizing Tacev fc Colel (|l993l , eq. 2.16), we find that 
the fraction of all 10* M© halos at z = 10 that are incor- 
porated into M > W^^Mq halos by z = 6 is 1.4 x 10"*. 
This fraction is much larger (by a factor of ~ 10'^) than 
the fraction of halos that have close and massive neighbors 
and are irradiated by a critical fiux. This is strongly sug- 
gestive that essentially all of these strongly irradiated halos 
are in the highly clustered tail of the spatial distribution, 
and will end up in massive halos, such as those hosting the 
SDSS quasars, by z = 6. For reference, we also note that 
the 1O*M0 halos contain a fraction 1.6 x 10~^ of the to- 
tal mass at z = 10, while IO'^^Mq halos at z = 6 have 
a mass-fraction of 2.8 x 10~^ at z — 6. This implies that 
^ (1.6 X 10"^) X (1.4 X 10"*)/(2.8 X lO"'') « 10% of the 
masses of the 10^^ Mq halos at 2 = 6 were already assem- 
bled into 10* Mq halos aX z — 10. In this picture, the z = 6 
supermassive black holes originally formed inside collapsing 
10* M© satellite halos of more massive star forming galaxies. 



4.5 Metal Enrichment of Tvir > 10'' K Halos 

While the above discussion shows that a non-negligible frac- 
tion of halos may be kept Il2-free by a strong LW back- 
ground, another possible effect that could invalidate the 
direct-collapse black hole formation scenario is the pres- 
ence of heavy elements in such halos. Indeed, most high- 
redshift DM halos with Tvir ~ lO'^K may be already enriched 
with at least trace amounts of metals and dust produced 
by pri or star-formation in their progenitors. lOmukai et al.l 
( j2008l ) recently studied the thermal and chemical evolu- 
tion of low-met allicity gas exposed to extremely strong UV 
radiation fields. They find that these metals and any ac- 
companying dust can catastrophically lower the gas tem- 
perature, producing an effective equation-of-state index of 
7 = d\np/d\n p <^ 1, even in the presence of a UV back- 
ground as strong as Jlw = 10^, which keeps H2 photodis- 
sociated. They suggest that gas fragmentation may thus be 
inevitable above a critical metallicity, whose value is be- 
tween Zci ~ 3 X 10~*Zq (in the absence of dust) and as 
low as ZcT ~ 5 X 10~^Zq (with a dust-to-gas mass ratio of 
about OMZ/Zq, also see Omukai 2000, Bromm et al. 2001, 
Schneider et al. 2003). When the metallicity exceeds these 
critical values, dense clusters of low-mass stars may form 
at the halo nucleus. While relatively massive stars in such a 
cluster could then rapidly coalesce into a single more mas- 
sive object (which may produce an intermediate-mass BH 
remnant with a mass up to M <^W'^ — 10^ M©), the forma- 
tion of much BHs as massive as M ^10^ - 10*^ M0 by direct 
collapse is likely to require - in addition to being exposed to 
a large Jcrit - that the halo gas is also essentially metal-free. 
Metal enrichment by the nearby star forming galaxy (that 
is responsible for the large local LW flux), is not necessarily 
a problem, since the expelled metals most efficiently mix in 
the outer edg e of the collapsing halo, while leaving the inner 



the outer edg e ot the collapsing halo, while leavmg the : 
core pristine l|Greif et al.ll2007l : ICen fc Riauelmdi2008l ) 
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5 CONCLUSIONS 

H2 molecules can be photo-dissociated by ultra-violet (UV) 
radiation, either directly (by photons with energies E > 14.7 
eV, if the molecules are exposed to ionizing radiation) or as 
a result of electronic excitation by Lyman- Werner (hereafter 
LW) photons with energies 11.2eV < E < En- In this paper 
we focused on this latter process, which operates even in 
gas that is self-shielded, and/or in gas that is shielded by a 
neutral intergalacti c medium (IGM) prio r to the completion 
of reionization (see iHaiman et al.lll997l ). from radiation at 
E > Eh- 

Photodissociation feedback possibly plays an impor- 
tant role in the formation of the supermassive black 
holes (SMBHs, Mbh ~ 10^ A/©) that existed at redshift 
z > 6, when the age of the universe was < 1 Gyr 
(e.g. lFanll2006l ). The reason for this is that this feedback 
mechanism can affect gas cooling and collapse into ha - 
los that have Tvir > lO" K (e.g. lOh fc HainianI l2002h . 
Without H2 molecules, gas inside these halos collapses 
nearly isothermally due to atomic cooling, with the tem- 
perature remaining as high as T ~ 10*K, which may 
strongly suppress the ability of gas to fr agment into stel- 
lar mass objects (e.g. lOh fc Haiman|[20o3 ). Instead of frag- 
menting, this gas could rapidly accrete onto a seed BH 
l|Volonteri fc Reej 120051 '). or collapse directly int o a very 
massive (Mbh = 10"^ - 10'^ Mp) b la ck hole ( B romm fc Loebl 



20031: 



2004 iBegeiinan et al.l I2OO6I : 



. .K oushiap pas et al.l . 

Spaans fc Silk.2006l : iLodato fc Nataraianll2006l . l2007l ). pos- 



sibly wit h an intermediate st ate in the form of a very massive 
star (see lOmukai et al]|2008l . for a more complete review). If 
these "direct-collapse black holes" indeed formed, then they 
would provide a head-start that could help explain the pres- 
ence of SMBHs with inferred masses of several xlO^M© by 

Possibly the most stringent requirement for the direct- 
collapse black hole model, is the intensity of the photodisso- 
ciating background, which must exceed J2i,lw ^ Jait ~ 10'' 
(Bromm & Loeb 2003, where J2i,lw denotes the intensity 
in the LW background in units of 10~^^ erg s~^ Hz~^ sr~^ 
cm~^, evaluated at the Lyman limit). Since the mean free 
path of Lyman- Werner photons is large, the LW-background 
is expected to b e nearly uniform, a nd likely significantly 
lower than Jcrit l|Omukai et al.l [20081 ). However, some vari- 
ation in the LW-flux that irradiates individual dark matter 
halos is expected. For example, a dark matter halo with a 
nearby star forming galaxy may 'see' a LW-flux that exceeds 
the background by orders of magnitude. 

In this paper, we have modeled this variation and we 
have computed the probability distribution function (PDF) 
of the LW-flux, Jlw, that irradiates a high-redshift dark 
matter halo aX z = 10. The main differences from previ- 
ous works that study the spatial variation of the UV back- 
ground are that: (i) we include (1) the non-linear clustering 
of sources, using the non-linear bias formalism developed by 
Iliev et al. (2003), which describes the clustering of closely 
separated halos particularly well; (2) Poisson fluctuations 
in the number of surrounding star forming galaxies; (3) the 
probable existence of a dispersion in LW luminosity for a 
given halo mass. In the local Universe the UV luminosity 
can vary by orders of magn itude for a given stellar mass 
(e.g lSchiminovich et al.ll2007l ). and a similar variation of the 



UV luminosity may be expected for a given halo mass; (ii) 
we specialize to compute the PDF of the flux as seen by 
high-redshift DM halos; (iii) we focus on the tail of the flux 
PDF, and (iv) we discuss the significance of this tail for 
rapid high-redshift SMBH growth. 

We found that > 99% of the dark matter halos at 2: = 10 
are illuminated by LW flux that is within a factor of 2 from 
the global background value (Fig. [Sjl . However, a tiny frac- 
tion (/crit ~ 10~* — 10~®) of dark matter halos with virial 
temperatures exceeding Tvir ~ lO** K may see a LW-flux 
that is boosted to the level Jcrit = 10'^ that is required by 
the direct collapse black hole model, due to nearby lumi- 
nous neighbors. We demonstrated that this result is insen- 
sitive both to our assumed model parameters (Fig. |6]), and 
to radiative transfer effects due to intervening atomic and 
molecular hydrogen gas (ij 14. ip . 

Depending on (i) the duration for which halos are pro- 
ducing UV luminosity (tium), (ii) the duration for which the 
irradiation is needed to avoid fragmentation (ttrag), and (iii) 
the spread in formation times of halos (tform), these rare halo 
pairs provide possible sites in which primordial gas clouds 
collapse directly into massive black holes (Mbh ~ 10''-^Mq, 
see i; l4.2|) . In 14.41 we showed that these black holes are likely 
to end up in massive (10^^ — 10^'' Mq) host halos (which also 
agrees with the mass of the host halos of the 2; ~ 6 quasars; 
see, e.g., the review by Haiman & Quataert 2004), in which 
they can grow to lif Mq within ~ 0.5 Gyr (the total time 
that elapsed between 2: = 10 and 2: = 6). It is interesting 
to stress that in this picture, the 2 = 6 supermassive black 
holes originally formed inside collapsing 10* satellite ha- 
los of more massive star forming galaxies. 

The fraction /cHt ~ 10"** -lO"'^ translates to ~ 10-10^ 
^10* Mq halos per comoving Gpc^ volume. The gas inside 
one of these halos needs to collapse directly into a massive 
black hole, and subsequently grow to a Mbh ~ 10^ Mq to 
explain observed space density of 2 = 6 luminous quasars. 
The actual number of supermassive black holes at 2 = 6 is 
higher by a factor of l/eqso, where eqso denotes the fraction 
of Hubble time that 2 = 6 supermassive black holes are 
observable as lum inous quasars. Assuming a quasar lifetime 
of ~ 50 Myr (e.g. iMartinilliooj '). we have eqso ~ 0.05, and 
the number density of 2 = 6 supermassive black holes is 
~ 20 cGpc"^, which we conclude is manageable with the 
direct-collapse black hole scenarios. 
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